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Abstract. We report spectro-imaging observations of the 
bright western ridge of the supernova remnant IC 443 ob- 
tained with the ISOCAM circular variable filter (CVF) on 
board the Infrared Space Observatory (ISO). This ridge 
corresponds to a location where the interaction between 
the blast wave of the supernova and ambient molecular 
gas is amongst the strongest. The CVF data show that 
the 5 to 14/Ltm spectrum is dominated by the pure rota- 
tional lines of molecular hydrogen (v = 0-0, S(2) to S(8) 
transitions). At all positions along the ridge, the H2 ro- 
tational lines are very strong with typical line fluxes of 
10~'*tolO~^ erg s~^cm~^sr~^. We compare the data to 
a new time-dependent shock model; the rotational line 
fluxes in IC 443 are reproduced within factors of 2 for 
evolutionary times between 1,000 and 2,000 years with a 
shock velocity of ~ 30kms^^ and a pre-shock density of 
- 10''cm-3. 

Key words: Infrared: spectra — ISM: supernova rem- 
nants — ISM: individual objects: IC443 — shock waves 



1. Introduction 

The supernova remnant IC 443 is a prime example of the 
interaction of a supernova blast wave with an ambient 
molecular cloud. On optical plates, IC443 appears as an 
incomplete shell of filaments (Fig. 1) with a total extent 
of about 20 arcmin, i.e. ^ 9 pc for an adopted distance of 
1500 pc. The shock generated by the supernova explosion, 
that occurred (4— 13) x 10"^ years ago, encountered nearby 
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molecular gas which is mainly found along a NW-SE di- 
rection across the face of the optical shell. IC 443 has been 
the subject of numerous studies from X-rays, visible, in- 
frared to radio wavelengths (e.g., Mufson et al., 1986 and 
references therein) . Studies of the interaction between the 
shock and the ambient molecular gas were done by observ- 
ing molecular hydrogen in the rotational-vibrational tran- 
sitions (Burton 1988, Burton et al. 1990- see Fig. 1 - and 
Richter et al. 1995a), in the pure rotational S(2) transition 
(Richter et al. 1995b), in other simple molecules such as 
CO and HCO+ (e.g., van Dishoeck et al. 1993 and refer- 
ences therein) and in atomic carbon (Keene at al. 1996). 
In this letter we report mapping results of the pure 
rotational lines of H2 using the ISOCAM CVF over the 
western ridge of IC 443, a position corresponding to clump 
G in the nomenclature of Huang et al. (1986). The obser- 
vations reveal the details of the structure and the physical 
conditions of the shocked molecular gas in IC 443 with a 
pixel field of view of 6" and at unprecedented sensitivity 
(mjy) . The observed H2 line fluxes are well predicted by a 
time-dependent shock model recently developed by Chieze 
et al. (1998) and Flower & Pineau des Forets (1999). 

2. Observations and data reduction 

The observations were done in 1998 February with the 
ISOCAM CVF (Cesarsky et al. 1996). A pixel size of 6" 
was used, yielding a total field of view of 3' x 3'. Scans of 
the long-wavelength CVF were obtained: the LW-CVF2 
from 13.53 to 9 fim and then the LW-CVFl from 9 to 
5.0 /im for a total of 115 wavelength steps; the resolving 
power is 40. Each wavelength was observed for 12.6 sec, 
i.e. six readouts at 2.1 sec, for a total observing time of 
some 30 minutes. 

ISOCAM was pointed towards a = 06h 13m 41.0s and 
S = 22° 33' 10.4" (coordinates B1950.0), a position corre- 
sponding to the center of the molecular clump G which is 
also a peak in the H2 emission (Fig. 1). 
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Fig. 1. The footprint of ISOCAM (box) depicted against the 
DSS image of the IC 443 supernova remnant together with the 
H2 l-OS(l) emission (contours) from Burton et al. (1990). 



The data reduction includes a new time dependent 
dark current correction (Biviano & Sauvage, in prepara- 
tion) as well as a new correction procedure for the tran- 
sients (Coulais & Abergel 1999). The photometric calibra- 
tion was done using the calibration files applicable to the 
current release of the ISOCAM off-line processing soft- 
ware (V7.0); the absolute calibration is conservatively es- 
timated to be on the order of 25%. 

The zodiacal background was subtracted adopting the 
spectrum obtained by Reach et al. (1996; but note that 
the published data were reduced again using the same 
algorithms and calibration files as applied to the current 
data) with a scaling factor of 0.7. 

The flux of the H2 lines was obtained by numerical in- 
tegration under the line profile, after subtracting a linear 
baseline defined by several points at either side of the line 
profile; the line fluxes thus obtained are independent of 
the assumed zodiacal background. We estimate the statis- 
tic uncertainty of the numerical integration by applying 
the same integration algorithm to regions of the spectrum 
devoid of any spectral emission; we consistently obtain an 
rms noise of some 2 10^^ erg s~^ cm~^ sr~^. 

3. Results 

The 5 to 14 /im spectrum of the molecular clump G in 
IC 443 is dominated by the series of the pure (w = — 0) 
rotational lines of molecular hydrogen from the S(2) to the 
S(8) transitions (Fig. 2). In particular, there is no indica- 
tion of any atomic fine structure line. At low level inten- 
sity, the set of dust emission bands from 6.2 to 11.3 /im is 
clearly present, including the 12.7 /xm band. Note that 
the relative strength of the dust bands is typical of that of 
the ISM (Boulanger 1998), suggesting that a possible con- 
tribution from the [Neil] line at 12.8 /zm is negligeable. 
This dust emission has comparable intensities over the en- 
tire ISOCAM field, i.e. a few MJysr"^, which is similar to 
the intensities measured in diffuse regions of the Ophiu- 
cus cloud. An image in one of the dust bands, e.g. 7.7 /im, 
does not show any structure correlated with the molecu- 
lar filament. Using the trend between the band intensity 



and the UV radiation field shown by Boulanger (1998), 
we find that the excitation of the dust bands is commen- 
surable with a radiation field G ~ 10 x Gq. The dust 
bands are probably unrelated to the supernova remnant 
and more likely mixed with the interstellar gas along the 
1.5 kpc line of sight toward IC443. Although very faint, 
the 7.7 and 6.2 /im dust bands contaminate the S(4) and 
S(6) H2 lines. We have fitted Lorentzians to the dust bands 
(Boulanger et al. 1998) and removed the resulting mean 
dust band spectrum from each individual CVF spectrum 
prior to performing the numerical integrations. 
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Fig. 2. fC443 mean CVF spectrum, i.e. integrated over 
the observed field, before (solid line) and after (dotted line) 
subtraction of the dust bands spectrum which is shown at 
the bottom of the figure at true scale but shifted down by 
4 10^^ erg s^^ cm^^ fim^^ sr^^. The rotational lines of molec- 
ular hydrogen are labeled. 



Figure 3 shows the total emission of the II2 lines be- 
tween 5 and 13.5 /xm, i.e. the sum of the S(2) to S(7) lines, 
together with the integrated line intensity of the S(2) and 
S(7) transitions (top panels). The II2 emission is found 
along a ridge of about 30" x 80" (0.25 pcx 0.65 pc) run- 
ning SW to NE, a structure which is comparable to that 
seen in CO or HCO+ (van Dishoeck et al. 1993, Tauber 
et al. 1994). The higher spatial resolution of the ISOCAM 
data clearly reveal a series of knots sitting on a plateau. 
The H2 knots are very bright with peak values of a few 
10"'^ erg s~^ cm~^ sr"-*^ (Table 1). The eastern side of the 
molecular ridge, facing the origin of the supernova explo- 
sion, appears sharper than the opposite side where weak 
emission is found extending westwards. 

4. Discussion 

Altogether there are about 130 pixels that show H2 emis- 
sion with intensities above the 10 a level, i.e. > 2 lO"** 
erg s^^ cm~^ sr^^ for all the six rotational transitions 
S(2) to S(7). This fact allows the construction of as many 
H2 excitation diagrams which plot the logarithm of the 
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Fig. 3. The distribution of the emission of the S(2) to S(7) H2 lines towards Clump G in IC443 (upper left panel). The next 
top panels show the emission in the S(2) and S(7) H2 lines. Contours are drawn at the 10%, 20%, etc. level. The corresponding 
peak strengths are 8.1 10"^, 7.9 lO""* and 2.2 10"^ erg s"^ cm~^ sr~^ from left to right. The middle panels show the ISOCAM 
CVF spectra towards the three emission peaks A, B, and C. The lower panels present the excitation diagrams after correcting 
for an extinction of A2.12 /im ~ 0.6 mag. Note that the spectra correspond to one pixel whereas the excitation diagrams have 
been built from a 3 x 3 pixel box car smoothed data. Each excitation diagram has been fitted with a two-component model 
involving 'warm' (dashed line) and 'hot' (dotted line) H2. The full lines represent the sum of both H2 components - see text 
and Tabic 1 for details. 



column density, corrected for statistical weight, in the up- 
per level of each H2 transition vs. the energy of that level, 
Eu. When calculating the column densities, we corrected 
the observed line fluxes for extinction, using a screen 
model and the extinction curve from Draine & Lee (1984). 
Peak A corresponds to a position studied by Moorhousc 
et al. (1991) and Richter et al. (1995a - their Position 3). 
Moorhouse et al. derived an extinction of A2.12 ^m = 0.6 
mag., whereas Richter et al. estimated a somewhat higher 
value (~ 1 mag.). We adopted the former value (see dered- 
dening factors in Table 1). The middle and bottom panels 
in Fig. 3 present the CVF spectra and the corresponding 
excitation diagrams of the three emission peaks labeled 
A, B, C. The statistical weights used in Fig. 3 include a 
factor of 3 for ortho-ll2 and 1 for para-H2. 

The excitation diagrams for Peaks A, B, and C show 
that a single excitation temperature does not reproduce 
the H2 lines observed in IC 443 and that emission from gas 
with a range of temperatures is required. The results of a 
simple LTE two-component H2 model are shown in Fig. 3: 
a 'warm' II2 component with an excitation temperature of 
~ 500 K and typical column densities in between 10^° and 
lO^^cm^^, and a 'hot' H2 component with Tex ^ 1200 K 
and Nh2 a few 10^^ cm^^ (Table 1). The parameters of the 
'warm' component are determined almost entirely by the 
intensities of the S(2) and S(3) lines and the 'hot' compo- 
nent dominates the H2 transitions S(4) to S(7). Although 
the evidence for a 'warm' component is very strong, the 
ISOCAM data only poorly constrain its properties be- 
cause of the lack of measurement of the S(l) and S(0) H2 
transitions. Furthermore, the uncertainty in the extinc- 
tion correction (especially for the S(3) line whose posi- 
tion coincides with the peak of the silicate 9.7 /xm band) 
introduces an additional uncertainty in the temperature 
determination. Using different extinction laws and adopt- 
ing values for A2.12 um between 0.5 and 1 mag., we derive 



typical uncertainties of ± 100 and ±250K for the warm 
and the hot components, respectively. 

ISO spectroscopy of regions where shocks dominate 
the excitation has revealed that the shocked gas has a 
range of temperatures from a few 100 K to several 1000 K: 
CepheusA (Wright et al. 1996), Orion (Rosenthal et al. 
1999) and bipolar outflows (Cabrit et al. 1999). The H2 
lines cannot be explained by a single shock model and 
combinations of J- and/or C-type shocks with different 
velocities and pre-shock densities have been invoked to 
account for the observed line fluxes. Similar conclusions 
have been reached for IC443. Based on an analysis of 
the [O i] 63 fim fine structure line and near-infrared H2 
lines. Burton et al. (1990) concluded that the infrared line 
emission of IC443 can only be modeled as a slow (10- 
20kms~^), partially dissociating J shock where the oxy- 
gen chemistry is suppressed, i.e. the cooling is dominated 
by [O l] emission and not by H2O cooling - see also Richter 
et al. (1995a, b). Far-infrared spectroscopy obtained with 
ISO on IC 443 confirms these conclusions and will be dis- 
cussed in a forthcoming paper. 

Following the interpretation of Wright et al. (1996) 
for the shocked H2 gas in CepheusA, the H2 lines in 
IC 443 can be fit by a combination of two C-shocks from 
the models of Kaufman & Ncufcld (1996). A good match 
to the data towards Peak A is obtained combining a first 
shock with a pre-shock density of 10'' cm^"^, a velocity of 
20kms^^ and a covering factor <i> of 0.85 with a second 
shock of lO^cm-^, 35kms-i and $ of 0.008 (sec Fig. 4). 
Such a steady-state model requires at least two C-shocks 
with a set of 3 free parameters and relatively high pre- 
shock densities for the high velocity component. 

Recently, Chieze et al. (1998) pointed out that the in- 
tensities of the ro-vibrational H2 lines are sensitive to the 
temporal evolution of a shock wave. In many astrophysi- 
cal situations, shock waves are unlikely to have reached 
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Table 1. Line fluxes and derived parameters from the model fits - see text 



H2 Transition 


S(2) 


S(3) 


S(4) 


S(5) 


S(6) 


S(7) 






Model fit 




Wavelength [fira] 


12.28 


9.66 


8.02 


6.91 


6.11 


5.51 


Warm Component 


Hot 


Component 


Dereddening factor^ 


1.16 


1.42 


1.13 


1.07 


1.08 


1.09 


T 


Nh, 




T 


Nh, 


Peak identification 


Observed line flux [10 


^■^ erg 


s^^ cm" 


-sr-1 


[K] 
657 


[lO^o cm 
2.2 


-'] 


[K] 
1288 


[10^0 cm-2] 


A 


3.7 


14.3 


7.7 


17.6 


4.4 


8.3 


0.22 


B 


5.1 


11.9 


5.3 


18.9 


5.1 


10.3 


330 


19.7 




1172 


0.48 


C 


2.4 


8.4 


4.5 


12.4 


3.3 


6.2 


446 


2.6 




1115 


0.37 



10' 



0.4 Aa 



adopting the extinction curve of Draine & Lee (1984) and A2.12 /im ~ 0.6 mag. 



steady-state which occurs at approximately 10^ yr. At 
times scales of a few 10'^ yr, the shocked gas may show both 
C- and J-type characteristics: within the C-shock, a J-type 
shock is established heating a small fraction of the gas to 
high temperatures (Flower & Pineau des Forets 1999). In 
the case of IC 443, the typical size of the molecular clump 
G is < 20", i.e. < 2 x 10^^ cm, but some of the molecular 
clumps have typical sizes of about a few arcsecs (Richter 
et al. 1995a). For a shock velocity of 20 — 30kms~^, the 
crossing time of the shock wave is thus < 2500 — 4000 yr 
indicating that the shock wave in clump G is not in steady- 
state. 




2 4 6 8 

E(v.J)/k (10^ K) 

Fig. 4. The excitation diagram for Peak A (filled circles) com- 
pared to the predictions of the time-dependent shock model of 
Chieze et al. (1998) with a pre-shock density of 10'* cm""', a 
shock velocity of 30kms~^ and four evolutionary times. The 
long-dashed curve labelled (KN96) presents the predictions of 
a model with two C-shocks from Kaufman & Neufeld (1996) - 
see text. 



Figure 4 shows the predictions of the time-dependent 
shock model for an observation along the direction of 
the shock propagation. The model results are given for 
four evolutionary times with the following parameters: 
pre-shock gas density uh = 10** cm^"^, shock velocity 
Vs = 30kms~^, magnetic field strength B == 100 /xG and 
ortho-to-para ratio of 3. The filling factor in this model 
is equal to 1. A smaller filling factor could be compen- 
sated by a larger line of sight path across the shocked 
H2 gas for a non face-on shock. The post-shock densities 
are 10"'' — 10® cm^^ comparable with the values derived 
by van Dishoeck et al. (1993). The agreement between 
the model predictions and the observations is best for in- 
termediate times, i.e. <2000yr. At earlier epochs when 
the J-shock dominates, the low excitation H2 lines are too 
weak. And after 5000 yr, when the C-shock steady-state 
is reached, the higher excitation H2 lines (above 10^ K) 
are much too weak. In between, the intensities of the H2 
rotational lines are predicted within factors of 2 and the 
coexistence of the 'hot' and 'warm' II2 components is well 
explained within a single model. Models with other param- 



eters (e.g., Uh = 3 X lO"^ cm 



35kms ) provide less 



good fits. The best fits are obtained for early evolution- 
ary times {^ 1000-2000 years) and densities of ^ 10^ cm""^ 
with shock velocities Vg = 30 — 40kms~^. Higher shock 
velocities will predict too large intensities for the high- 
excitation II2 lines. 
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The agreement between these time-dependent model 
predictions and the data of IC 443 is very encouraging in 
view of some of the simplistic underlying assumptions. In 
particular, the assumed geometry (plane parallel) is over- 
simplified and does not describe the molecular filament 
which is seen edge-on and consists of numerous small 
clumps. Clearly a more thorough study should be done 
to explore the entire parameter space of the model and 
compare the predictions with additional data available on 
IC443. 

In the model, the S(2) to S(7) lines account for about 
70% of the luminosity in all the II2 lines. At peak A, the 
measured S(2) to S(7) flux is ^ 4.6 x 10^^^ erg s^^ cm^^ 
(~ 0.4 Lq) and, according to the model, the H2 lines alone 
would thus carry ~ 0.6Lq. Towards clump G, the S(3) and 
S(2) H2 lines account for almost the entire IRAS 12 ^m- 
band emission. The mean value of these lines in that band 
is 5.4MJysr~^ comparable to the IRAS peak value, i.e. 
6MJysr-i (e.g., Oliva et al. 1999). Similarly, the IRAS 25 
/im-band could also be due to II2 line emission. The model 
predictions (Fig. 4) for the S(0) and S(l) line fluxes (at an 
evolutionary time of 2000 years) are 1.2 10~^ and 7.7 lO^"' 
erg s~^ cm~^ sr~^, respectively. Taking into account the 
20% transmission at 17.03 /xm of the IRAS 25 /im band, 
the strong S(l) line would thus contribute ^ 4MJysr^^ at 
25 iJ,m., which is comparable to the measured 25 /im IRAS 
flux at Peak A (- 4.5MJysr-\ e.g. Oliva et al. 1999). 
These results strongly suggest that the excitation of the 
gas in IC 443 is entirely collisional. 

Finally, Oliva et al. (1999) found towards the optical 
filaments of IC 443 , which trace the low density atomic 
gas, that most of the 12 and 25 /im IRAS fluxes is ac- 
counted for by ionized line emission (mainly [Neil] and 
[Fell]). Our results show that this conclusion cannot be 
generalised towards the molecular hydrogen ring (Fig. 1) 
where the dense molecular gas essentially cools via the 
H2 lines in the near- and mid-infrared and via the [O i] 
emission line in the far-infrared. 

Acknowledgments: Michael Burton is kindly thanked 
for providing his II2 map of IC 443. 
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